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Abstract. The mean motion of turbulent patterns detected by a two-dimensional (2D) beam 
emission spectroscopy (BES) diagnostic on the Mega Amp Spherical Tokamak (MAST) is de- 
termined using a cross-correlation time delay (CCTD) method. Statistical reliability of the 
method is studied by means of synthetic data analysis. The experimental measurements on 
MAST indicate that the apparent mean poloidal motion of the turbulent density patterns in the 
lab frame arises because the longest correlation direction of the patterns (parallel to the local 
background magnetic fields) is not parallel to the direction of the fastest mean plasma flows 
(usually toroidal when strong neutral beam injection is present). The experimental measure- 
ments are consistent with the mean motion of plasma being toroidal. The sum of all other 
contributions (mean poloidal plasma flow, phase velocity of the density patterns in the plasma 
frame, non-linear effects, etc.) to the apparent mean poloidal velocity of the density patterns 
is found to be negligible. These results hold in all investigated L-mode, H-mode and inter- 
nal transport barrier (ITB) discharges. The one exception is a high-poloidal-beta (the ratio of 
the plasma pressure to the poloidal magnetic field energy density) discharge, where a large 
magnetic island exists. In this case BES detects very little motion. This effect is currently 
theoretically unexplained. 

(Some figures in this article are in colour only in the electronic version) 

PACS: 28.52.-s, 52.55.Fa, 52.70.Kz, 52.30.-q, 52.35.Ra, 52.35.Kt 
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1 Introduction 

It is now widely accepted that turbulent transport in magnetically confined fusion plasmas can 
exceed the irreducible level of neoclassical transport by an order of magnitude or more [UJ. 
However, both theoretical and experimental works of the past two decades [12H15D suggest that 



sheared E x B flows can moderate such anomalous transport and hence improve the perfor- 
mance of magnetically confined fusion plasmas. 

With the aim of characterizing the microscale plasma turbulence and searching for corre- 
lations between it and the background plasma characteristics, a two-dimensional (8 radial x 



4 poloidal channels) beam emission spectroscopy (2D BES) system 111611 has been installed on 
the Mega Amp Spherical Tokamak (MAST). It is able to measure density fluctuations at scales 
above the ion Larmor radius pi, viz., kj_pi < 1, where k± is the wavenumber perpendicular to 
the magnetic field. The 2D BES view plane lies on a radial-poloidal plane at a fixed toroidal 
location. Following the detected turbulent density patterns on this view plane allows one to 
determine their mean velocity in the radial and poloidal directions. Typically, there are no sig- 
nificant mean plasma flows in the radial direction in a tokamak, whereas considerable apparent 
poloidal motion is detected by the 2D BES system. 

In this paper, we show that this apparent poloidal motion is primarily due to fact that the 
direction of the longest correlation of the turbulent density patterns is not parallel to that of 
the dominant mean plasma flows. The BES measurements are shown to be consistent with a 
dominantly toroidal mean flow; the poloidal flows are of the order of the diamagnetic veloci- 
ties. These results are obtained using the cross-correlation time delay (CCTD) method, which 
is a frequently used statistical technique to determine the apparent velocity of density patterns 



fll7l ll8H. We also investigate the method itself thoroughly to determine the statistical uncertain- 
ties of the technique. This is done by generating synthetic 2D BES data with random Gaussian 
density patterns calculated on a graphical processing unit (GPU) card using CUDA (Compute 
Unified Device Architecture) programming. 

The paper is organized as follows. In section |2l we explain what is measured directly by 
the 2D BES system, and how the apparent velocity of turbulent density patterns can be in- 
ferred from this data. We also show what physical effects contribute to the apparent velocity 
calculated by the CCTD method. In section [3l we present the experimental results with the 
aim of identifying the main cause of apparent motion of density patterns measured by the 2D 
BES system. Our conclusions are presented in section |4] For the readers who are interested 
in the statistical technique employed in this paper to determine the velocity of the density pat- 
terns, the cross-correlation time delay (CCTD) method and its statistical reliability are studied 
in Appendix B using synthetically generated 2D BES data (described in Appendix A[ ). 



2 What is measured by the 2D BES system 



The 2D BES system on MAST utilizes an avalanche photodiode (APD) array camera [|19ll 
with 8 radial and 4 poloidal channels, which have an active area of 1.6 x 1.6 mm 2 each. It 
measures the Doppler- shifted D a emission from the collisionally excited neutral beam atoms 
(deuterium) with a temporal resolution of 0.5 psec. The optical system is designed so that the 
2D BES system can scan radially along the neutral beam, whose 1/e half-width is 8 cm, while 
the optical focal point follows the axis of the beam. The nominal location of the BES system, 
i.e., where the optical line-of-sight (LoS) is best aligned with the local magnetic field, is at 
major radius R = 1.2 m. At this location, a magnification factor of 8.7 at the axis of the beam 
results in each channel observing an area of 1.5 x 1.5 cm 2 with 2 cm separation between the 



Motion of turbulent density patterns detected by the BES system on MAST 4 

centres of adjacent channels. The angle between the LoS of the 2D BES system and the neutral 
beam with the injection energy of 60 — 70keV results in a Doppler shift of the D a emission 
approximately 3 nm to the red from the background D a . The background D a can be removed 
with a suitable optical filter, and so only the detected D a emission by the 2D BES system 
comes from the neutral beam, hence the locality for measurement to the beam. Aligning the 
LoS parallel to the local magnetic field at the intersection of the LoS and the neutral heating 
beam helps minimize the degradation of the spatial resolution. A more detailed description of 
the 2D BES system on MAST with possible sources of some losses of spatial locality can be 



found elsewhere IU6U20I1 . 



2.1 Plasma density fluctuations 

The measured intensity of the D a beam emission is directly related to the background plasma 
density because the latter is the cause of the excitation of the neutral beam atoms. The beam 
atoms can be excited by electrons, ions and impurities, but for the injection energy greater than 



40 keV, the electron contribution can be ignored [12111 . The fluctuating part of the plasma (ion) 
density 5n can be determined according to 

n pi 

where n is the mean plasma density, 51 and / denote the fluctuating and mean parts of the 
photon intensity, respectively, and /3 is calculated based on the ADAS (Atomic Data and Anal- 
ysis Structure) database [|22Jj. f3 is a weak function of the background plasma density ranging 
approximately from 1/3 to 1/2. 

Thus, the 2D BES system on MAST directly measures fluctuations of plasma density in the 
poloidal-radial plane at a fixed toroidal location. The spatial resolution of the system is reduced 
by smearing due to the effects of field-line curvature, observation geometry, finite lifetime of 
the excited neutral-beam atoms, and the attenuation and divergence of the beam. These effects 
must all be taken into account in the calculation of the point spread functions (PSFs) of the 
detectors comprising the 2D BES system. A detailed calculation shows that ~ 2 — 3 cm. radial 
resolution and ~ 1 — 5cm poloidal resolution are achievable, depending somewhat on the radial 



viewing locations Il23h 



2.2 Velocity of density patterns 

From the time-dependent 2D measurement of density fluctuations, one can infer th e ap parent 
velocity of the density patterns. This has been the subject of much attention Ill7l.l24l-l29ll in the 
hope that this velocity can be related in a more or less straightforward way to the actual plasma 
flows. We will first explore how the mean pattern velocity can be determined and then discuss 
the interpretation of this quantity. 

2.2.1 The CCTD method 

The CCTD (cross -correlation time delay) method has been widely used to determine the appar- 
ent velocities of turbulent density patterns detected by BES systems, and it is well described in 



ll 1 Vh and II 1811 . Here, a brief summary of the method is provided. The normalized fluctuating 



intensity of the photons, I = 51 / 1, measured by a 2D BES system is a function of the radial x, 



Motion of turbulent density patterns detected by the BES system on MAST 





1.0 


(a) A xv : 


6 


[(b) 


""''; 






'h\\ 

: / ' / V w ■ 


5 


D = 11.4 ±0.1 [A:/wA] / 


/ - 





0.5 


1 i \\\\ 






: 


H 




: <<\ / / / , :\v\' j 


4 




- 


5 




" <, x \\ /'/'/' 'il f - 


" 




: 


o" 


0.0 


\\\ //' \ \ ^ ' 1 1 ' 


3. 3 




: 


II 










- 


3 




1 V*\- /'• WW 7 - 7 / 


2 




: 


(J 


-0.5 


XT 1 ' 1 \ x y 






: 






A 1 l Ay-O.OM 


1 




j 






•V / 2.0c»i 






; 






4.0cm — — 






; 




-1.0 


6.0cm — - — - 







- 




-20 -10 10 20 


12 3 4 5 


6 






x [jisec] 




Aj [c/w] 





Figure 1. (a) Cross-correlation functions calculated using equation Q for Ay = 0.0 cm (black solid 
line), 2.0cm (red dash line), 4.0cm (blue dash dot line) and 6.0cm (green dash dot dot line). T c p c eak (Ay) 
is the position of maximum of the cross-correlation function, (b) Position of maximum r^ ak (Ay) and 
a lineal - fit. The measured velocity is 11.4 ± 0.1 km/s. 



vertical (poloidal) y and time t coordinates: I = I (x, y, t). The cross-correlation function of 
this fluctuating signal is defined as 



C(Ax,Ay,r) 



/ (x,y,t) I (x + Ax,y + Ay,t + r) 



P {x,y,t)) (P (x + Ax,y + Ay,t + r] 



(2) 



where Ax and Ay are the radial and vertical (poloidal) channel separation distances, respec- 
tively, r is the time lag, and (■) denotes time average defined in |Appendix B[ The apparent 



poloidal velocity v^ ES of the density patterns detected by the 2D BES system can be deter- 



mined from the time la£ 
a given Ay and Ax = C 



r. 



peak 



at which the cross -correlation function reaches its maximum for 



If a straight line is fitted to the experimentally measured T^ ak (Ay), 



'peak 

the inverse of its slope is the velocity v^ ES . Although any two poloidally separated channels 



are sufficient to determine v BES , using just two channels is insufficient to estimate the uncer- 
tainties in the linear fit. Thus, in this paper, all four available poloidal channels are used to 
determine these quantities. This assumes that the mean velocity does not change over the time 
the density patterns take to move past the four poloidal channels and that the lifetime of these 
patterns is sufficiently long, so the same patterns are observed by all four channels. 

Figure \T\ shows an example of this procedure. This example is based on a synthetic data set 
consisting of Gaussian-shaped random "eddies" moving with the poloidal velocity of 10.0km/ s, 



which are then used to produce artificial 2D BES data (see Appendix A for the description of 
the synthetic data). With the four available poloidal channels, cross-correlation functions are 
calculated using equation © and shown in Figure [TJa); r^ ak is plotted as a function of Ay in 
Figure QIb). The inverse of the slope of a fitted straight line is the velocity v BES . Note the 



*We concentrate on the apparent mean 'poloidal' motion of the density patterns. Thus, the information about 
the radial correlations of the 2D BES data is not used in this paper. 
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Figure 2. Cartoon illustrating how the mean toroidal plasma flow (U z ) induces an apparent mean 
poloidal motion. An elongated density pattern (shaded oval) along the magnetic field line (green dash 
dot) is advected by the toroidal flow (blue arrow). Because the longest correlation direction of the 
density pattern is not in the toroidal direction, the apparent mean poloidal flow (green arrow) arises. The 



apparent velocity is —U z tan a + U y 



-U z tan a, where a is the local magnetic pitch angle. 



slight discrepancy between the actual and CCTD-determined velocities. The origin and size of 



this discrepancy are discussed in Appendix B 



2.2.2 Physical meaning of the velocity determined by the CCTD method 

Using the described CCTD method, the 2D BES system on MAST is expected to be able to 
determin e v B ES as has been done previously on TFTR Ill7ll and DIII-D 112911 using their BES 
systems IBOL Bill . However, as McKee et al. II32L 13311 pointed out, one must distinguish be- 
tween the poloidal velocity measured by 2D BES system (v BES ) and the actual velocity of the 
poloidal plasma flow (U y ). 

The mean plasma flow can be decomposed into toroidal (U z ) and poloidal (U y ) components. 
For typical tokamak plasmas where strong neutral beams are injected, \U Z \ 3> \U y \ is satisfied 
as any mean poloidal flows are strongly damped [B4l 13511 . leaving U y of the order of the dia- 
magnetic velocity ~ p*v th , where p* = Pi/a, a is the tokamak minor radius, and v th is the ion 
thermal velocity. Note that U z can be on the order of v t h for the neutral-beam-heated plasmas. 
Thus, U y can be ignored compared to U z , except possibly in regions with strong pressure gra- 
dients. 

As the 2D BES system on MAST observes the density patterns advected by U z , there will 
be an apparent motion of the patterns in the poloidal direction, as shown in Figure [2] This 
effect is analogous to the apparent up-down motion of helical strips of a 'rotating barber-pole' 
(cf. [36;]). The magnitude of this apparent velocity can be readily calculated via elementary 
geometry: namely, we expect the BES system to "see", to lowest order in p*, 



,BES 



—U z tana, 



(3) 
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where a is the pitch angle of the local magnetic field line. 

Equation © is experimentally verifiable because all three physical quantities are readily 
obtained by separate diagnostics: v BES from the 2D BES system, U z from the charge exchange 

J I 

recombination spectroscopy (CXRS) system [37], and a either from EF IT equilibrium recon- 



struction [1381] or directly from the Motional Stark Effect (MSE) system [39L [40Q on MAST. 
Although the CXRS system measures the toroidal flow of the C 6+ ions, the difference between 
the velocity of the C 6+ ions and the bulk plasma ions, D + , is predicted to be on the order of 



p* in a strongly beam-heated plasma 114 ill . In section [3j equation © will be experimentally 
verified for various types of discharges. Agreement will indicate consistency of the experiment 
with the assumptions behind equation ©. Such agreement will indeed be obtained, except in 
one intriguing case. 

Let us now consider what are the assumptions necessary for equation ([3]) to hold by analysing 
how the estimated v BES depends on actual physical quantities associated with plasma flows and 
fluctuations in a tokamak. The cross -correlation function © of the normalized fluctuating pho- 
ton intensity I can, in view of equation CD, be considered proportional to the cross -correlation 
function of the relative ion density fluctuation Sn/n (by definition, (5n) = 0). Therefore, the 
CCTD-determined velocity of the density patterns can be related to the actual physical quanti- 
ties in a tokamak by invoking the ion continuity equation. Splitting also the ion velocity into 
mean and fluctuating parts, u = U + Su, (Sit) = 0, we have 



dn dSn 



V • (nU + nSu + SnU + 8n5uj = 0. (4) 



dt dt 
Averaging this equation and subtracting the averaged equation from ©, we obtain 

dSn 



Dt 



-V • ( nSu + SnU + SnSu - (SnSu)) . (5) 



We will now order various terms in this equation in terms of the small parameter p* = Pi/a. 

Assuming that the spatial scale of all mean quantities is ~ 0(a) while the spatial scale of 
all fluctuating quantities is ~ 0(p*a), and also 5n/n ~ 5u/v t h ~ p*, we get 

= -17 -V V-Su 

at n n 

—5u)- — (v-U + U-V\nn)+0(pl), (6) 

where we have dropped all terms ~ C(p#) and smaller. The first two terms on the right-hand- 
side are ~ 0(v t h/a) and the following three terms are ~ 0(p*v t h/a). Note that we have not 
yet made any assumptions about the nature of the mean flow U (beyond it being large-scale) or 
about time scale of the fluctuations. 

In fact, the p* ordering, which is the standard gyrokinetic ordering [14211 . can take us further: 
it is possible to show that compressibility effects are order p*, i.e., V • Su ~ 0(p*), and that the 



mean flow to lowest order is purely toroidal 1134113 511 : U = U z z + U\, where z is the toroidal 



direction (locally) and U\ ~ @(p*) including all poloidal flowi] and first-order corrections 
to U z (radial flows, associated with particle fluxes, are, in fact, even smaller). Coupled with 
the fact that mean quantities have no toroidal variation in a tokamak, this means that that the 



^Note that the poloidal velocity U y of the bulk plasma ions has been measured with the CXRS system to be 



only a few km/ s on MAST 14311 . which is consistent with U y ~ 0(p*). Such measurements are, however, not 
routinely available for MAST, and one of the goals for this study is to confirm that U y is indeed small. 
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fifth term on the right-hand-side of equation © is also ~ O(pl), while the first term can be 
expressed as 

T -t „Sn TT d 5n r - T „Sn 

U-V— = U z - + U 1 -V— 

n oz n n 

= -l/,^ + -r^b - V— + Ui ■ V— , (7) 

b z ay n b z n n 

where b = (0, b y , b z ) is the unit vector in the direction of the magnetic field in a local orthog- 
onal Cartesian system (x: radial, y: poloidal and z: toroidal), and we have used the identity 
b ■ V = byd/dy + b z d/dz. Making a further assumption, again standard in gyrokinetics, that 
the parallel spatial scale of the fluctuating quantities is ~ 0(a), we conclude that the second 
term in the second line of equation © is 0(p*). 
Finally, combining equation © with ©, we find 

d 5n TX d 5n 5n 

— — + U eff - = 7 —, (8) 

at n ay n n 

where U e fj = —U z b y /b z = —U z tan a is the dominant apparent velocity of the density patterns 
(a is the local pitch angle of the magnetic field line). The term containing U e ff is the only 0(p®) 
term in equation ([8]). The 0(p*) and higher terms have been assembled in the right-hand-side: 
by definition, 7 is such that 



5n 


U Z 7 ^Sn 


7 — = 


: —±b-V — 


n 


b z n 



[5n r 



U x - V— 

n 

-V -Su-Su-Vlan-V ■ i—5u\ + O(p'i). (9) 

This contains, in order of terms, the effects associated with 

(1) parallel variations of the fluctuations, 

(2) mean poloidal flows of bulk plasma ions, 

(3) compressibility of the fluctuations, 

(4) linear response to mean density gradient (drift waves), 

(5) nonlinear effects (turbulence), 

and a slew of higher-order effects of varying degree of obscurity. 

Thus, the right-hand- side of equation © contains all the nontrivial physics of waves and 
turbulence in the plasma. The apparent velocity of the density patterns detected by the 2D BES 
system will not be influenced by these effects to dominant order — if the orderings assumed 
above are correct. What it does contain is the poloidal signature [/ e // of the dominant toroidal 
rotation of the plasma — the 'rotating barber-pole' effect discussed at the beginning of this 
section. Indeed, if equation © holds and its right-hand- side is small, then, to lowest order, the 
density patterns just drift in the y-direction (poloidal) with the velocity U e ff, so the maximum 
of the cross-correlation function © will be achieved at r = Ay/U e ff. Hence equation © for 
the BES-measured velocity. 

If we are able to confirm equation © experimentally, this means that the theoretical consid- 
erations employed above are consistent with the experiment. This is important because most of 
the theories of tokamak turbulence rely on such considerations. Note that there are no separate 
diagnostics capable of measuring individually all the 0(p*) terms in equation ®. Therefore, 
the only conclusion one can formally draw from equation © holding is that the sum of these 
terms is small. 
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3 Experimental results 

In this section, we apply the CCTD method to 2D BES data from MAST discharges to de- 
termine the apparent mean poloidal motion {v EES ) of the ion density patterns. Then, v EES is 
compared with the 'rotating barber-pole' velocity (U z tan a) where the toroidal plasma velocity 



fci 



U z is obtained from the CXRS system 07h and the local magnetic pitch angle a either from 



EFIT equilibrium reconstruction Q38J or the MSE system []39l 14011 . 

The 2D BES data are first bandpass-filtered from 20.0 to 100.0 kHz to reduce the noise 
level. The low-pass filter removes the high-frequency noise component from the photon shot 
noise and electronic noise, while the high-pass filter reduces the contribution to the signal from 
low-frequency, coherent MHD (magnetohydrodynamic) modes. The apparent mean poloidal 
velocity of the density patterns v BES is determined from average correlation functions calcu- 
lated over 25 time intervals of 40 fxsec duration, resulting in total 1 msec averaging. Second- 
order polynomial fitting is applied to interpolate the correlation function on times shorter than 



the sampling time of 0.5 usee as described in Appendix B.l Finally, five consecutive values 
of v BES obtained in this manner are averaged, so the total averaging time is 5 msec which is 
the effective time resolution of Vy ES . Using these five values of v EES , the time average of 
various errors defined in equations (fl~8l>-(l20l in Appendix B.2| are also computed. Statistical 



reliability of the CCTD method is investigated in Appendix B by using the synthetic 2D BES 



data (Appendix A). 



We present measurements of v BES from four different discharges: shot #27278 (L-mode), 
shot #27276 (H-mode), shot #27269 (ITB) and shot #27385 (high-poloidal-beta). All four dis- 
charges had double-null diverted (DND) magnetic configurations and co-current NBI (neutral 
beam injection). In all of these discharges, the 2D BES system viewed at nominal major ra- 
dial position of R = 1.2 m corresponding to normalized minor radii r/a = 0.2 — 0.3 for L- 
and H-modes, and r/a = 0.3-0.4 for ITB and high-poloidal-beta discharges. The evolution 
of key parameters for these discharges is shown in Figure |3] The evolution of plasma cur- 
rent, line-integrated electron density and poloidal beta characterize the overall behaviour of 
plasmas, while the non-zero S-beam voltage corresponds to times when the 2D BES system 
obtains localized density fluctuation. The D a intensity trace is used to identify when the H- 
mode discharge (shot #27276) goes into its H-mode: namely, at t = 0.21 — 0.28 sec. Note that 
the ITB discharge (shot #27269) starts developing a strong temperature gradient at ~ 0.2 sec 
and the peak ion (C 6+ from the CXRS) temperature keeps increasing until the NBI cuts off at 
0.3 sec. The viewing position of the 2D BES system is in the middle of the strong temperature 
gradient region for this discharge. 

3.1 L-mode (shot #27278), H-mode (shot #27276) and ITB (shot #27269) 
discharges: Vy ES « — £4 tan a 

Time evolution of (a) cross-power of the fluctuating magnetic field signal from two toroidally 
separated outboard Mirnov coils, (b) cross-power and (c) temporal cross -correlation of den- 
sity fluctuations from two poloidally separated BES channels (two mid-channels separated by 
2 cm) located at R = 1.21 m are shown in Figures |4] (L-mode), [5] (H-mode) and [6] (ITB dis- 
charge). Here, a cross-power is defined as the Fourier transform (in the time domain) of the 
cross-correlation function © with finite channel separation. The (minus) apparent mean 
poloidal velocity (—Vy ES , circles) determined by the CCTD method and the 'rotating barber- 
pole' velocity (U z tana, red solid lines) are also shown in panels (d) at R = 1.13 m and (e) at 
R = 1.21 m for these three discharges. The error bars represent the mean error (Svf it ) of the 
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Figure 3. Evolution of (a) plasma current (b) line-integrated electron density (c) poloidal beta (d) NBI 
(S-beam) injection energy and (e) edge D a intensity of L-mode (shot #27278, black solid), H-mode 
(shot #27276, red dash), ITB (shot #27269, green dash dot) and high-poloidal-beta (shot #27385, blue 
dash dot dot) discharges. 



least-squares fit, as discusses in Appendix B.2 



Despite the fact that these three discharges belong to three very different classes, there are 
common features in the apparent mean poloidal velocity: 

(1) Vy ES is not reliable (i.e., has large error bars) when strong MHD activity is present. The 
cross-power spectrograms from BES show clear signatures of MHD modes with many har- 
monics, which hamper filtering the BES signal over the frequency domain. The temporal cross- 
correlations also show that these MHD modes have much longer correlation times (> 0.3msec) 
than the turbulent density patterns. The effects of MHD (global) modes on the CCTD method 
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Time |sec| 

Figure 4. The evolution of shot #27278 (L-mode) showing (a) cross-power spectrogram of the fluctu- 
ating magnetic field signal from two toroidally separated outboard Mimov coils, (b) cross-power spec- 
trogram and (c) cross-correlation of the density fluctuations between two poloidally separated channels 
(two mid-channels separated by 2 cm) from BES at R = 1.21 m. The time evolution of the (minus) 
apparent mean poloidal velocity (—Vy, circles) from BES and the 'rotating barber-pole' velocity 
(C/ztana, red solid line) from the CXRS at (d) R = 1.13 m and (e) R = 1.21 m. BES signals are 
bandpass-filtered from 20.0 - 100.0 kHz for (c)-(e). 



are investigated in Appendix B.5 where it is found that such activity can increase not only the 



absolute values of the bias errors but also the linear fitting errors on Vy ES '. Thus, comparisons 
between v BES and U z tan a are difficult to make during the periods where the MHD activity is 
strong. 

(2) During the periods of weak MHD activity, i.e., 0.11-0.15 sec for the L- and H-mode dis- 
charges, and 0.16-0.22 sec for the ITB discharge, it is clear that the apparent mean poloidal 
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Figure 5. Same as Figure|4]for shot #27276 (H-mode). 



velocity of turbulent density patterns is dominated by the 'rotating barber-pole' velocity, i.e., 
equation © holds, and the sum of all the terms of the order of p* or higher in equation © is 
indeed small. 

Note that the H-mode discharge (shot #27276) goes into its H-mode at ~ 0.21 sec (thus, 
v BES = —U z tan a is only true before the L-H transition, strictly speaking), which can be seen 
from the D a intensity trace in Figure |3]or from the BES cross-power spectrogram in Figure |5] 
the turbulence level drops at the start of the H-mode. Any changes of Vy ES during the L-H tran- 
sition cannot be discussed, because the CCTD method with the current data analysis scheme is 
not reliable at this time due to strong MHD activity. 
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Figure 6. Same as Figure [4]for shot #27269 (ITB discharge). 



3.2 



High-poloidal-beta discharge (shot #27385): v% ES ^ 



■UA&na 



Shot #27385 has a relatively higher poloidal beta (the ratio of the plasma pressure to the 
poloidal magnetic field energy density) than the three discharges discussed in section I3TI (see 
Figure 3j. Thus, it is more susceptible to tearing modes (i.e., formation of magnetic islands) 
11441. 14511. The cross-power spectrogram between the two toroidally separated outboard Mirnov 
coils displayed in Figure |7£a) shows a m/n = 3/2 tearing mode on the g=1.5 flux surface start- 
ing at ~ 0.11 sec; its frequency increases from < 10 kHz to ~ 25 kHz at ~ 0.19 sec. Then, a 
m/n = 2/1 mode (fundamental frequency < 10 kHz) develops and locks to the wall resulting 
in complete braking of the toroidal rotation of plasmas at ~ 0.25 sec. Here, m and n denote 
the poloidal and toroidal mode numbers, respectively, and q for the safety factor. 

The 2/1 mode is not expected to be seen on the BES signal as it is bandpass-filtered from 
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Figure 7. Same as Figure |4] for shot #27385 (high-poloidal-beta discharge). Note that neither (b) 
the cross-power spectrogram nor (c) the temporal cross-correlation of the BES signal show any MHD 
activity. 



20.0 — 100.0 kHz, and no trace of the 3/2 mode is visible in the BES signao Consequently, 
the v BES determined by the CCTD method does not contain large error bars during the whole 
discharge. 

The time evolution of —v BES and U z tan a in Figure [2d)-(e) at two different radial loca- 
tions shows that the two velocities do not agree each other at all during the period when the 3/2 
and 2/1 modes are present. What we find, remarkably, is that while the plasma continues to 
rotate toroidally (as attested by the CXRS data), there is virtually no detectable corresponding 



^Because the mode flattens the mean density profile within the island, shaking of flux surfaces does not induce 
density fluctuations in the BES signal 
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motion of the density patterns. In fact, they seem to exhibit a weak rotation in the opposite 
direction to the expected rotating barber-pole effect. Formally, this means that the plasma ef- 
fects in the right-hand-side of equation ([8]) are not small and are able to cancel almost exactly 
the toroidal rotation, i.e., an effective velocity of the density patterns develops in the plasma 
frame that to lowest order is equal to minus the rotation velocity. We do not currently have 
a theoretical explanation for this effect. There is very little apparent difference between the 
turbulent density patterns in this discharge compared to others, except somewhat longer radial 
correlation lengths. 

4 Conclusions 

We have analysed 2D BES data from different types of discharges on MAST to determine the 
apparent mean poloidal velocities of the ion-scale density patterns using the cross-correlation 
time delay method. The dominant cause of the apparent poloidal motion of the density patterns 
is experimentally identified to be due to the fact that field aligned patterns are advected by the 
background, dominantly toroidal, plasma rotational flow, i.e., the 'rotating barber-pole' effect 
dominates the apparent mean motion of the density patterns in the lab frame. This conclusion 
holds for the L-, H-mode and ITB discharges we have investigated. An exception to this rule 
is found to be the investigated high-poloidal-beta discharge, where a large magnetic island is 
present, and the apparent velocity of the density patterns is very small, despite strong toroidal 
rotation. Identifying the causes of this effect by investigating the behaviour of the turbulent 
density patterns quantitatively is left for future work. 
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Appendix A Synthetic 2D BES data 
Appendix A.l Gaussian eddies in space and time 

It is necessary to know the true mean velocity of the density patterns to investigate statistical 
reliability of the cross-correlation time delay (CCTD) method described in section [272] Specif- 
ically, we investigate how reliable the CCTD method is for different magnitudes of mean ve- 
locities and correlation times of the density patterns. We must also evaluate the effect of global 
modes (i.e., MHD modes) and temporally varying poloidal velocities on the CCTD method. 
For this purpose, we numerically generate artificial density patterns, random both in space and 
time, then produce synthetic BES data using the point- spread-functions (PSFs) of the 2D BES 



system on MAST (as described in Appendix A. 2 ) and compare the inferred flow velocity with 
the true flow velocity. 
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We follow a similar approach to the one suggested by Zoletnik et al. II46I1 . Let the density 
patterns be described by Gaussian structures both in space and time, namely, 



N 



Sn (x, y,t) = ^ Sn oi ex P 



cos 



(x-x 0i ) 2 [y + Vy (t) (t - tot) -y i] 2 (t-t 0i ) 2 ' 



2vr 



2A 2 2 A 2 2r 2 

x y life 

[y + v y (t) (t - t 0i ) -y 0i ] 

\y 



where x, y and t denote radial, poloidal and time coordinates, respectively. These numerically 
generated density patterns are referred to as "eddies" in this paper. Here N is the total number 
of eddies and the subscript i denotes the i th eddy in the simulation; 5n 0i , x 0i , y 0i and t 0i are 
the maximum amplitude and central locations in the x, y and t coordinates of the i th eddy, 
respectively; \ x , \ y and ruf e are the widths of our Gaussian eddies in the x, y and t directions; 
Tuf e is the lifetime (or the correlation time) of the eddies in the moving frame; v y (t) is the 
apparent advection velocity of the eddies in the poloidal direction. Although it is possible to 
introduce a finite radial velocity shear by making v y a function of x, the effect of such shearing 
rates on the CCTD method is not investigated in this paper, so we will only consider v y that 
are independent of x. The cos term in the y (poloidal) direction is introduced to model wave- 



like-structured eddies in the poloidal direction as observed in tokamaks 114711 . Note that the 
envelope (i.e., the exp term) and the wave structure (i.e., the cos term) of 5n(x, y, t) have the 
same advection velocity v y (t). The central locations of eddies, x , y and t , are selected 
from uniformly distributed random numbers, whereas their amplitudes 5n are selected from 
normally distributed random numbers whose standard deviation is one|j 

The spatial domain of the simulation is 25 cm and 20 cm with the mesh size of 0.5 cm 
in radial (x) and poloidal (y) directions, respectively. The time duration of the simulation is 
20 msec with a 0.5 jjsec time step so as to have the same Nyquist frequency as the real 2D 
BES data from MAST. The widths X x and X y are set so that the full width at half maximum 
(FWHM) in the radial direction and the wavelength in the poloidal direction are ~ 8 cm (i.e., 
\ x = 3.53cm) and ~ 20cm (i.e., \ y = 20.0cm), respectively, which are similar to the measured 
correlation lengths with the 2D BES system on MASTjjThe eddy lifetime in the moving frame 
(rufe) is set to 15 usee. However, some of the data sets in this paper have different values of 
Tuf e , so the effect of Tuf e on the CCTD method can be investigated. 

The total number of eddies is N = 20000. If the eddies are too sparse in the simulation 
domain, then we may not achieve steady statistical results, while overly dense eddies may cause 
an effective widening of the specified spatial (A x and X y ) and temporal (rufe) correlations as 
many eddies can merge into one larger eddy. Thus, we introduce another control parameter, the 
spatio-temporal filling factor (F), defined as 

p _ jy . ( * x *y ] . ( Jac \ ,jj, 

\ total simulation area / V total simulation time/ ' 



§It is worth mentioning that there is another scheme of generating such eddies numerically, proposed by 



Jakubowski et al. 14811 . They generated the time series of fluctuating density (8ni) using the inverse Fourier 
transform of a broadband Gaussian amplitude distribution in frequency space. Then, a second signal (^2) was 
generated by imposing the desired time-delay fluctuation on the 5n\ such that 5ri2 was a time-delayed version of 
8n\. This method does not include spatial information for the signals. 

^Note that Smith et al. fl49fl also reported that poloidal correlation lengths of the density patterns are ~ 20 cm 
using their 2D BES system on NSTX. 



(10) 
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where r ac is the autocorrelation time calculated as 

„ _ T K f e (Xy/Vy) 



2 _i_{\ /„, \ 2 



Tu 



life 



(Xy/Vy 



for the generated eddies defined by equation (HOI ). All of our synthetic data was generated so 
astoF~ 0(1). 

The testing of the CCTD method will involve exploiting what happens if v y (t) has a mean 
and a temporally varying components. Thus, we generate a temporal structure of v y : at each x, 



Vy(t) = (Vy) + Sv y (t) 

= (v y ) + v y (t) * exp 



t 2 



T fluc 



m(2irf fluc t) (13) 



sin 



where (v y ) and 5v y are the mean and temporally varying velocities, respectively, Tf iuc and ff tuc 
are the lifetime and frequency of Sv y (t), respectively, and v y (t) is generated from normally dis- 
tributed random numbers. The RMS (root-mean- square) value of 5v y (t) denoted as 5Vy MS will 
be varied as well as (v y ) to investigate the effects of these quantities on the CCTD method. Tfi uc 
and ff Luc allow one to introduce structured temporally varying velocities, while the randomness 
is kept by v y . As one of the causes for the temporal variation of the poloidal velocity is believed 
to be the existence of geodesic acoustic modes (GAMs]U Q5O0, we choose Tfi uc = 500 usee and 
ffiuc = 10 kHz to mimic the GAM features detected by Langmuir probes on MAST pill . 

The simulations have been run on a NVIDIA® GeForce GTS 250 GPU card using CUDA 
programming, which increases the computational speed owing to the highly parallelizable 
structure of equation (fTOl . 

Appendix A.2 Synthetic 2D BES data 

We generate the i th (1 to 8) radial and j th (1 to 4) poloidal channel of the synthetic BES data 
Pi (t) byusing the calculated point- spread-functions (PSFs) of the actual 2D BES system on 



MAST ll23n and Sn (x, y, t) from equation (flOl) with an additional random noise. Furthermore, 
a large-scale (in space) coherent (in time) oscillation is included to imitate a global MHD mode. 
Namely, P J (t) is defined as 

I lj (t) = I% C + SPi (t) + I% HD (t) +!%{t), (14) 



where I % f] c is the DC value - a typical value of 0.8 V is used for all channels [16]. The rest of 
the terms are as follows. 

51 % i (t) is the fluctuating part of the signal generated from the Gaussian eddies, 5n (x, y, t), 
given by equation (fTOl) and convolved with the PSFs of the 2D BES system : 




5P 3 (t) = 5I KMb / / 5n (x, y, t) V lJ (x, y) dxdy, (15) 



where P JJ (x, y) is the PSF of the i th and j th channel of the 2D BES system, normalized so that 
RMS value of SP j (t) is 5I RMS . This value is set so that the ratio of 5I RMS to I% c is 0.05. An 
example of the PSFs for the 32 channels of the 2D BES system on MAST is shown in Figure [9] 



II We do not investigate whether the CCTD method is able to detect such a temporally structured 8v y (t) (or 
GAMs) in this paper, rather we investigate how the existence of these structures affects the CCTD-determined 
mean velocity. 
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Figure 8. Autopower spectra of synthetic 2D BES data for various (v y ). Note that the spectrum for 
(v y ) = 10.0 km/ s (green dash dot line) has finite I^ IHD (i.e., temporal oscillations due to global 
modes) in equation (fT4l) at 15 kHz, with fluctuation level of 5 % of the DC level. For other cases, 
I mild = 0- AH the spectra are generated using a high-pass filter with the frequency cutoff at 5 kHz. 



(white contour lines in the top-left panel). 

1-mhd (0 m °dels an MHD (global) mode. We assume that the spatial scale of the MHD 
modes is larger than the BES domain in the poloidal direction, so P* 



mhd (0 does not vary 

in the poloidal direction. The model MHD signal is generated in a way similar to temporal 
behaviour of v y (t) using equation (fl~3T) . except that the mean value of I MHD (t) is zero and 
Tfi uc = 250 usee. The frequency of the mode /mhd and its RMS value, denoted I^ H S D , 
will be varied in various tests. The value of Tn U c here is representative of MHD burst-like 
fishbone instabilities Il52ll or chirping modes 115311 in tokamaks, for which the spectrum has a 
finite bandwidth. 

(t) represents the noise in the signal. As the noise of the 2D BES system on MAST 



1 N 



is dominated by the photon noise IU9I1 . 1% (t) is generated using normally distributed random 



numbers. Its RMS level is set such that the signal-to-noise ratio (SNR) is 300, which is typical 
of the 2D BES system on MAST H]. 

Figure [8] shows examples of autopower spectra of the synthetic 2D BES data for (v y ) = 
2.0, 5.0, 10.0 and 40.0 km/s. The autopower spectrum is calculated as \FT{I l i (t)}\ 2 where 
FT {■} is the Fourier transform in the time domain. Increasing the value of (v y ) has two effects: 
Doppler shift and broadening of the spectra, as expected. Note that in Figure [81 the data for 
(v y ) = 10.0 km/s contains the finite I^mj with } M hd = 15 kHz and I^hdI^dc = °- 05 ' 



while I mhd — f° r other cases. 



Figure|9]shows several time snapshots of artificial Gaussian eddies (equation (fTOl) ) in the left 



Motion of turbulent density patterns detected by the BES system on MAST 



19 



Pseudorandom eddies 




21) synthe tic BES data 

>*t > • 'i ■ 






At=20us 



At=30^is 







1 .00 



1.05 t.10 1.15 
.v Jm] (Major Radius) 



1.20 



1.00 1,05 1.10 1.15 
x [in] (Major Radius) 



1.20 



Figure 9. Left column: four time snapshots of Gaussian eddies, 5n (x, y, t) given by equation (flOl) . 
Right column: the corresponding normalized synthetic 2D BES data given by equation (fl4t w ithout the 
DC component. White lines in the top left panel show the 1/e contour lines of the PSFs J23I1 . and the 
white asterisks show the optical focal points of the 32 channels of the 2D BES system. 



column and the corresponding synthetic 2D BES data in the right column (with DC component 
removed from equation (O). The eddies are moving upward with (v y ) = 5.0 km/ s. The top 
left panel in this figure also shows the 1/e contour lines of the PSFs for the 32 channels [23]. 
Snapshots for the synthetic 2D BES data are generated with the bandpass frequency filtering 
from 10 to 70 kHz to suppress the noise. As the synthetic 2D BES data have only 32 spatial 
points, spatial interpolation is performed using parametric cubic convolution technique Il54ll . 



Appendix B Assessment of the CCTD method 

In this section, errors involved in determining the mean velocity of the density patterns by the 
CCTD method are examined using the synthetic 2D BES data generated according to the pro- 



cedure explained in Appendix A The velocity measured via the correlation function (equation 
©) is denoted v^ ES and compared with the prescribed value (v y ) that appears in equation 



131) . i.e., the mean poloidal velocity of the synthetic data. Appendix B.l provides detailed 



description of the CCTD method used in this paper, then four types of error are identified for 
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the quantitative comparisons. These errors are evaluated in Appendix B.3 and Appendix B.4 
for different values of (v y ) and the eddy correlation time Tuf e . Subsequent sections are devoted 
to investigating how the existence of global (MHD) modes and temporally varying poloidal 
velocity affect the errors. 

Appendix B.l Description of the CCTD method 

As defined by equation ©, cross-correlation functions are calculated as time averages of the 
data. For a 20msec- long synthetic data set containing N to tai = 40, 000 data points with the sam- 
pling time At sam = 0.5 usee, we want to determine v BES with a time resolution t res = 1 msec. 
First, a cross-correlation function © is calculated on a sub-time window of the synthetic 2D 
BES data containing N c points, where N c < t res /At sam . Then, such cross -correlation func- 
tions are averaged over N avg consecutive sub-time windows where N avg = (t res /At sam )/N c 
so that an averaged cross-correlation function is obtained at every t res . In this paper, we use 
N c = 80, so N avg = 25. 

Denoting /(£) and g(t) the time series over a sub-time window from two poloidally sepa- 
rated synthetic 2D BES channels, the cross-correlation function © for this sub-time window 
is: 

Nc-l 

it 5Z fik^tsam) g{(k + r)At sam ) 



L'sub y^-^sam) 



k=0 



N C 



N c -l N c - 



(16) 



J2f(kAt sam ) Y,9 2 {{k + r)At 

\ k=0 



samj 
k=0 



for any integer r with \r\ < Nc — 1. Finally, by averaging C su b for N avg consecutive sub-time 
windows we obtain the smoothed averaged cross-correlation function C(rAt sam ) from 1 msec- 
long data points. 

The CCTD method has a serious limitation due to the fact that the sampling time At sam 
is finite. In order to calculate v BES using only two poloidally separated channels, a line is 
fitted through two points on a (Ay, Tp^ ak ) plane as shown in Figure Qlb). The first point is 
located at (Ay,r^ ak ) = (0,0) by definition, and the second point at (Ay,r At sam ). Then, 
possible values of v BES are restricted to Ay/ (rAt sam ) where r is an integer. For the 2D BES 
system on MAST, using two adjacent poloidal channels (Ay = 2.0 cm) with a sampling time 
At sam = 0.5 jj,sec, the possible values of v BES are limited to 40.0, 20.0, 13.3, . . . km/ s for 
r = 1, 2, 3, . . .. Such a limitation may be mitigated by using four poloidally separated channels. 
However, using four channels is not always possible if the channels that are farthest apart are not 
correlated. To resolve this issue, we use a second-order polynomial fit on the cross -correlation 
function C(rAt sam ) to locate its global maximum: if r peak is the point where the discrete cross- 
correlation function C(rAt sam ) is maximum, we use the three values of C(r At sam ) at r = r peak , 
r P eak — 1 an d r peak + 1 to fit a second-order polynomial. The "true" maximum is found from 
this fit. We denote the time delay at which this maximum is reached by r^ ak . 

Appendix B.2 Definition of errors 

For a given set of 20 msec-long synthetic 2D BES data, we calculate v BES with the time reso 



lution of 1 msec (Appendix B.l ). Furthermore, we do this at three different radial location 



**As described in | Appendix A| v y (t) are identical at all radial locations. One column in the middle and two 
columns from the edges of the 2D BES channels are used. 
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so that the average of v BES , denoted (v BES ^j, can be calculated using 60 values of v BES . To 
make quantitative comparisons between (v BES ) and (v y ) defined in equation (fT3l) . we define 
four types of error. 

The normalized bias error 

&bias = : — v (1 /) 

( v v) 

is a quantitative measurement of the systematic discrepancy between the measured and the true 

value. The normalized random error 



&rand \~i opcx i \*-°) 

quantifies the degree of fluctuation in the measured v BES with respect to (v BES ). This value 
may depend on the MHD contribution in equation (fl4l) and the temporally varying poloidal 
velocity 5v y (t) in equation (PT3T) . 

Furthermore, as linear fitting is done to determine v BES (see Figure [D, two other types of 
error are present. The slope of a linear fit can be denoted as v BES ± Svf it where 5vj it is a degree 
of the uncertainty of the least-square fl© Then, the normalized mean of Svfu is 

&f it = (Svfit) (19) 

" mean / brc\ ' ^ ' 

\\ v y )\ 
and the normalized random error in 5vf it is 



-fit 



yiifofit - (5v fit )f) 



\{v BES )\ 

These two uncertainties together provide an estimation of how well a linear line is fitted to given 
data points. For example, if the assumption that r^/e is long enough so that all four poloidally 
separated channels observe the same eddies is not satisfied, then o^an becomes large. On the 
other hand, if this assumption is occasionally satisfied, then (J r ^ nd exhibits such events because 
Sv fa will then be small compared to its average. Note that error bars of the CCTD -determined 
apparent velocities in Figures |4]-|7] show (5vf it ). 

In the following sections, these four types of error will be evaluated for various values of 
(v y ) and r life , and various ranges of Imh S d , Jmhd and SvP*' 



u 



Appendix B.3 Measuring mean velocity 



To investigate the reliability of the CCTD method described in Appendix B.l for estimating 



v 



BES 



we generate a number of synthetic 2D BES data sets with various values of (v y ) while 
keeping all the other parameters in equations (fTOl i. (fT3l ). (fT"4l) and (fT5l) constant. In real exper- 
iments, there is almost always some temporal variation of v y , thus the RMS value of Sv y in 
equation (PT3T) is set to 5 % of (v y ) in this subsection. The synthetic 2D BES data are frequency- 
filtered to suppress the noise before the cross-correlation functions are calculated. Figure \\0\ 
shows examples of (a) v y (£) generated according to equation <TT3T > with (v y ) = 10.0 km/ s and 
(b) the original (black) and frequency-filtered (red) autopower spectra of a generated synthetic 
signal. Here, the noise cut-off level is set to be the 5 times the averaged autopower level above 
900 kHz (green dashed line). 



ft In FigureQ] we plotted r"t k as a function of Ay and determined Vy ES as the inverse of the slope of a fitted 



line. Operationally, we actually plot Ay as a function of T™ fe so the slope of a fitted line is the Vy ES . 
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Figure 10. (a) Poloidal velocity v y (t) generated using equation (PT3T ) with (v y ) = 10.0 km/s. (b) 
Autopower spectra of the original (black) and frequency-filtered (red) synthetic BES signals. The green 
horizontal dashed line shows the noise cut-off level, defined to be 5 times the averaged autopower level 
above 900 kHz, and vertical blue dash-dotted lines indicate the low- and high-frequency cutoffs. 
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Figure 11. The four types of error defined in equations dT7)-(|20} calculated for values of (v y ) ranging 
from 1 to 100 km/s. 



Figure [TT] shows a bi 



Vrand, o f JL n and a 



fit 



'rand defined in | Appendix B.2| and calculated for 
values of (v y ) ranging from 1 to 100 km/s. The basic conclusions that can be made based on 
these results are as follows: 

(1) For (v y ) < 5.0 km/s, the CCTD method is not reliable. This is due to the fact that eddies 
do not live long enough to be detected by all the poloidally separated channels. Indeed, it was 
a priori clear that (v y ) < Ay/ruf e could not be measured. This translates to (v y ) < 4.0 km/s 
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Figure 12. Relative discrepancy between the propagation time r prop = Ay/ (v y 

(black) or T^™ k (red) at which the cross-correlation function or its envelope reaches their peaks for (a) 

(v y ) = 5.0 km/s and (b) 20.0 km/s. 



for Ay = 6.0 cm and Tnf e = 15.0 \isec, so our results are consistent with this simple criterion. 

(2) The CCTD method usually overestimates (v y ) (i.e., atias > 0). This can be explained by 
the effective channel separation distance (Ay) being in fact slightly less than 2.0 cm because of 
the overlapping of the PSFs, as shown in Figure [9] 

(3) The limitation of the CCTD method due to the finite At sam is successfully overcome by fit- 



ting a second order polynomial to the cross -correlation function, as explained in Appendix B. 1 



Appendix B.4 Effect of the eddy lifetime 



As explained in section 12". 2 .2[ the CCTD method for determining (v y ) is based on the idea that 
the peak of the cross -correlation function occurs at r^ ak = r prop , where r prop = Ay/ (v y ) is 
the propagation time of the fluctuating density patterns between detectors poloidally separated 
by the distance Ay. However, r^ ak will not coincide with r prop if the lifetime Tuf e of the 
fluctuations is not long compared to r prop . The failure of the method for (v y ) < 5.0 km/s 
illustrated in Figure \TT\ is an example of what happens when T prop is too large. Here, we 
investigate the effect of ruf e on T™ k quantitatively, via a systematic Tuf e scan of the synthetic 
BES data. 

Two values (v y ) = 5.0 and 20.0 km/s are chosen for this study. For (v y ) = 5.0 km/s, 
Tprop = 4.0, 8.0 and 12.0 usee with Ay = 2.0, 4.0 and 6.0 cm, respectively; for (v y ) = 
20.0 km/s, they are 1.0, 2.0 and 3.0 usee. The peak time T™ k is found using the polynomial 



fitting method described in |Appendix B.l[ and (r prop — rfg ofe ) /t ptop as a function of r^j e is 



plotted for three different values of Ay in Figure [[2l It shows that r^ eak underestimates the 
true Tprop for small values of Tuf e , leading to an overestimation of the (v y ), consistent with the 
results shown in Figure QT] It is encouraging, however, that even relatively low velocities of 
just a few km/s can be determined by the CCTD method with reasonable accuracy (~ 20%). 
It is also possible to consider the global maximum of the envelope of the cross-correlation 
function. We use Hilbert transform to determine the time delay rT^ k at which the envelope 
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Figure 13. Four types of error (a) as functions of the RMS levels of a global mode iff^D relative to 
that of turbulence signal 5I RMS ; the frequency is fixed at /mhd = 15.0 kHz; (b) as functions of the 
global mode frequency /mhd at fixed I R ^ S D /5I RMS = 5.0. Note that auas in (a) is scaled down by a 
factor of 10, and some points are missing in (b) because they are out of the plot range. 



of the cross-correlation function is maximum 11171) . The comparison between r^ k and T prop is 



shown in Figure Q21 It is clear that T^™ k has a much stronger dependence on ruf e than r^ ak , so 
this measure will not be used to estimate (v y ) in this paper. We note, however, that the strong 
dependence of rf™ fe on the eddies' lifetime Tuf e and of r^ ak on their propagation time T prop 
may provide a way to measure correlation times in the plasma frame. Such an investigation is 
currently being pursued and will be reported elsewhere. 



Appendix B.5 Effect of coherent MHD modes 

Many experimental 2D BES data sets on MAST exhibit strong MHD (global mode) activity 
in addition to the small-scale turbulence. Removing such global modes in the frequency do- 
main is not straightforward as they can have multiple harmonics extending into higher frequen- 
cies. While they could be filtered out relatively easily in the wavenumber domain, constructing 
wavenumber spectra with a very limited number of spatial data points is difficult. Thus, it is 
useful to investigate how the presence of such modes affects the quality of our measurement of 
(v y ). In this section, this is done by using synthetic BES data sets with different RMS levels 
Im A hd an d frequencies /mhd of the global oscillations (the I^hd term m equation (fl4l). 

The four errors (a bias , a rand , fr£? ean and a^ nd ) are calculated for various ratio of Imhd 
to the RMS value of 5P j (t) (i.e., SI RMS in equation (Tl3T)). These errors are plotted in Fig- 
ure USa) for the Iff^f) scan. Here, the frequency of the global mode /mhd = 15.0 kHz 
and (v y ) = 10.0 km/ s. It is clear that if the power level of the mode is larger than that of 
the turbulence signal, then the CCTD method produces large bias errors (r bias . To examine 
how the frequency of a global mode affects the errors, /mhd is varied with a fixed value of 
Imhd/3I RMS = 5-0. The results of this scan are shown in Figure [T3T b). It shows that au as can 
be either positive or negative with different values of /mhd meaning that global modes in real 
experimental data can cause both over- and under-estimation of the true (v y ). 
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Figure 14. Cross-correlation functions of the random eddies only (red dash), the global mode only 
(blue dash dot) and the eddies with the global mode (black solid) with (a) /mhd = 15.0 and (b) 
/mhd = 40.0 kHz. Green arrows indicate the position of T^ ak , which does not coincide with the 
maximum of the cross-correlation function of the eddies only (red dash). 



Figure[T4]shows how different frequencies Jmhd can cause such an over- or under-estimation 
of the (v y ). Two identical sets of synthetic BES data with (v y ) = 10.0 km/ s are generated, one 
with and another without a global mode at (a) Jmhd = 15.0 kHz and (b) /mhd = 40.0 kHz, 
with Imh£)/5I rms = 5.0. Without the global modes, the cross-correlation functions with 
Ay = 6.0cm (red dashes in Figured!]) have the expected value r^ ak rs 6.0 /isec for both cases. 
In contrast, the presence of the global mode in the synthetic BES data shifts r^ eak towards (a) 
smaller time-lag (over-estimation) or (b) larger time-lag (under-estimation). 

We conclude that a global (MHD) mode with ij^f, > SI RMS affects the structure of 
the cross-correlation functions (both the shape and the position of r^ ak ) rendering the CCTD 
method unreliable. 



Appendix B.6 Effect of temporally varying poloidal velocity 

No physical quantities are absolutely quiet in real experiments, thus it is necessary to investigate 
how the RMS level 5v RMS of the temporal variation of the poloidal velocity (see equation (fT3l) ) 
influences the measurement of (v y ). 

Figure [15] shows how finite Sv RMS / (v y ) (with (v y ) = 10.0 km/s) affect the four errors 



defined in Appendix B.2 It appears that auas saturates at around 50% for the scenarios we have 
investigated, while other three errors increase without showing any sign of saturation. Thus, 
the CCTD method to measure (v y ) is subject to a non-negligible bias error (up to ~ 50%) if 
the RMS level of temporal variation of the poloidal velocity is greater than a half of the mean 
poloidal velocity. 
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